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En este trabajo se presentan resultados del estudio de un 
humedal de flujo subsuperficial horizontal con Cortadeira 
selloana construido en el año 2010 para tratar el efluente 
cloacal generado en una dependencia de la Universidad 
Nacional de Córdoba, Argentina. La construcción del mismo 
se realizó con el objetivo de analizar la factibilidad de tratar 
los efluentes cloacales de una vivienda tipo, por lo que se 
construyó a escala real con un área superficial de 10 mY, 
aproximadamente. Se calculó el valor de la constante cinética 
de remoción de materia orgánica y se diseñó un software 
de simulación que permite evaluar distintas situaciones 
hipotéticas considerando al sistema como un reactor flujo 
pistón (REP). También se realizó un ensayo con trazador 
conservativo para determinar la distribución del tiempo 
de residencia (DTR) y el comportamiento hidrodinámico 
del flujo. Como resultado, se pudo detectar que el sistema 
en estudio no se comporta según el modelo hidrodinámico 
supuesto para su diseño, sino que se asemeja a un sistema 
conformado por dos reactores flujo pistón con dispersión 
en paralelo. Se obtuvo una constante de reacción para la 
eliminación de materia orgánica, considerando un modelo 
de REP y una cinética de remoción de primer orden de 
k = 0.115 d*. La simulación permitió conocer a priori la 
concentración de materia orgánica del efluente del humedal 
para situaciones de interés. Se destaca la importancia de 
realizar estudios en sistemas de una escala mayor a pruebas 
pilotos, tal como se realizó en el presente trabajo (en el cual 
el sistema se hizo a escala real), y monitorearlos en el tiempo, 
ya que las raíces de la vegetación, acumulación de sólidos, 
asentamiento de la grava y depresiones en el suelo por el 
peso de la celda son factores que hacen que el sistema se 
modifique y cambie su morfología y por lo tanto cambie su 
capacidad depurativa. 


Palabras clave: humedal construido, Cortadeira selloana, 
comportamiento hidráulico, constante cinética. 


Abstract 


Alasino, N., Nadal, F., Primo, L., Bollo, P., €: Larrosa, N. (May-June, 
2015). Kinetic and Hydraulic Behavior of a Full-Scale Constructed 
Wetland. Water Technology and Sciences (in Spanish), 6(3), 
93-104, 


This work presents the results from the study of horizontal 
subsurface flow in a constructed wetland with Cortadeira Selloana. 
The wetland was built in 2010 to treat waste generated by a facility 
belonging to the National University of Cordoba, Argentina. It 
was constructed to analyze the feasibility of treating wastewater 
generated by housing, and was therefore built at full-scale with 
an area of approximately 10 m?. The kinetic constant for removing 
organic matter was calculated and simulation software was designed 
to evaluate different hypothetical situations considering a piston 
flow reactor system. A conservative tracer test was also performed to 
determine the distribution of residence time and the hydrodynamic 
behavior of the flow. As a result, it was determined that the system 
does not behave according to the hydrodynamic design model but 
rather as a system composed of two piston flow reactors with parallel 
dispersion. A reaction constant was obtained for the removal of 
organic matter based on a piston flow reactor model and a first-order 
removal kinetics of k= 0.115 d*. The simulation identified a priori 
the concentration of organic matter in the effluent from the wetland 
for the situations of interest. The importance of studies using scales 
that are larger than pilot tests is highlighted (as was done in the 
present study using a full-scale system), as well as monitoring them 
over time, given that the vegetation roots, accumulation of solids, 
settling of gravel and soil depressions from the weight of the cell 
change the system and its morphology and, therefore, its purification 


capacity. 


Keywords: Constructed wetland, Cortadeira selloana, hydraulic 
behavior, kinetic constant. 


Recibido: 07/08/2013 
Aceptado: 12/02/2015 


Tecnología y Ciencias del Agua, vol. VI, núm. 3, mayo-junio de 2015, pp. 93-104 


gua, vol. VI, núm. 3, mayo-junio de 2015 


Ciencias del A 


« vy80/0u72], 


Alasino et al., Comportamiento cinético e hidráulico de un humedal construido a escala real 


Introducción 


Los humedales artificiales son sistemas de fito- 
depuración de aguas residuales. El sistema con- 
siste en el desarrollo de un cultivo de macrófitas 
de hojas flotantes o enraizadas sobre un lecho 
de grava impermeabilizado. La acción de las 
macrófitas hace posible una serie de complejas 
interacciones físicas, químicas y biológicas, a 
través de las cuales el agua residual afluente es 
depurada progresiva y lentamente (Delgadillo, 
Camacho, Pérez, € Andrade, 2010). 

El tratamiento de efluentes por humedales 
puede ser descrito como un tratamiento biológi- 
co y microbiológico, donde la materia orgánica 
de las aguas residuales es metabolizada por 
las bacterias del sistema. El oxígeno necesario 
para el metabolismo aerobio de estos microor- 
ganismos es aportado por la fotosíntesis de las 
plantas a través de sus raíces y rizomas. La sedi- 
mentación (en regiones de flujo lento), adsorción 
en las partículas de suelo y en los sistemas de 
raíces y otras partes de plantas, fotooxidación, 
captación directa de la planta e intercambio 
catiónico también son importantes mecanismos 
de remoción de contaminantes (Woodard and 
Curran, Inc., 2011). 

Los humedales han cobrado relevancia en la 
ingeniería sanitaria debido a que son sistemas 
que permiten una muy flexible localización 
de sitio, optimización del tamaño según la 
carga contaminante y tratamiento de una gran 
cantidad de aguas residuales. Constituyen una 
solución de bajo costo ampliamente utilizada en 
el mundo para el tratamiento secundario y/o 
terciario de diferentes tipos de aguas residuales 
(Vymazal, 2009). 

Los humedales construidos para el trata- 
miento de agua residual se han clasificado de 
manera tradicional en dos tipologías, aten- 
diendo a si la circulación del agua es de tipo 
subterránea o superficial. 

En los humedales de flujo superficial (HFS) 
el agua está expuesta directamente a la atmós- 
fera, con una profundidad máxima de 0.4 m y 
circula de modo preferente a través de los tallos 
de las plantas. 


En los humedales de flujo subsuperficial 
(HFSS), la circulación del agua es subterránea, 
a través de un medio granular con una profun- 
didad de alrededor de 0.6 m, y en contacto con 
los rizomas y raíces de las plantas que crecen 
gracias a los espacios de la porosidad de la 
grava. Los humedales de flujo subsuperficial 
se clasifican según el sentido de circulación 
del agua en horizontales (HFSSH) o verticales 
(HESSV). 

Los humedales, durante el transcurso del 
tiempo, pueden cambiar las condiciones bajo 
las cuales fueron diseñados y construidos, 
cambiando su morfología y, con esto, su con- 
ductividad hidráulica y porosidad, incidiendo 
de forma notable en los resultados esperados. 

En este trabajo se presentan los resultados 
de estudios realizados en un humedal construi- 
do —con Cortadeira selloana en el mes de junio 
de 2010— para tratar los efluentes cloacales 
generados en una dependencia de la Facultad 
de Ciencias Exactas, Físicas y Naturales de la 
Universidad Nacional de Córdoba, Argentina, y 
de una vivienda familiar situada de manera con- 
tigua al mismo (-31* 26' 21.29”, -64? 11' 28.47”). 

Para el diseño del humedal y cálculo de sus 
dimensiones se utilizó la metodología de Reed, 
Miledlebrooks y Crites (1995), que considera al 
humedal como un RFP con cinética de degra- 
dación de materia orgánica de primer orden 
(ecuación (1)): 


DBO), efluente _ pe 
DBO, afluente (1) 


donde DBO, efluente y DBO, afluente son las 
concentración de demanda bioquímica de oxí- 
geno (DBO,) en el efluente y afluente, respecti- 
vamente; k, la constante de velocidad cinética, 
y tes el tiempo de retención hidráulica. Para el 
cálculo de 1 se aplicó la ecuación (2): 
Mal 0) 
v 
donde L es el largo; W, el ancho; h, la profundi- 
dad del humedal; e, la porosidad del lecho, y v 
es el caudal volumétrico del efluente. 
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Se consideró para su diseño un caudal (v) de 
7001/d (consumo de agua promedio de las dos 
instalaciones al momento de la construcción). El 
resto de los valores se tomaron según recomen- 
daciones bibliográficas. Se consideró una DBO, 
afluente de 240 mgO,/l; DBO, efluente de 30 
mgoO, /1 (valor exigido por la DIPAS (1999) para 
descarga directa en embalse); e de 0.38; kde 0.73 
d! (Lara-Borrero, 1999); h de 0.5 m, y relación 
L/W=3.El área resultante de este humedal fue 
de unos 10 m?, con 2.1 m de ancho y 5.1 m de 
largo. El caudal de diseño y volumen de la celda 
arrojó un tiempo de residencia hidráulico teóri- 
co de 2.9 días (Nadal, 2010; Nadal et al., 2012). 

Durante los ocho primeros meses desde 
su construcción (desde julio de 2010 hasta fe- 
brero de 2011), el humedal trató agua residual 
doméstica con un caudal promedio de 700 1/d 
proveniente de las instalaciones sanitarias de la 
dependencia de la universidad y la vivienda. A 
partir de marzo de 2011, la vivienda fue desha- 
bitada y el caudal de ingreso al humedal dismi- 
nuyó a unos 222 1/d promedio. En la figura 1 se 
muestran fotografías del humedal construido. 


Materiales y métodos 


Durante los años 2011 y 2012 se llevaron a cabo 
diversos estudios que permitieron evaluar el 


funcionamiento del sistema: comportamiento 
hidrodinámico, evaluación cinética y simula- 
ciones de situaciones hipotéticas mediante un 
software diseñado para tal fin. 


Ensayo de trazador para determinación de u 
experimental y evaluación hidrodinámica 


La evaluación hidrodinámica es una herramien- 
ta útil para conocer el comportamiento que 
describe el sistema en la práctica y así ajustarlo 
al modelo de reactor o combinación de reactores 
que mejor se aproxime. Para ello es preciso rea- 
lizar ensayos con trazadores no conservativos, 
con el fin de obtener la curva de distribución de 
tiempos de residencia (DTR) (Hincapié-Marín, 
2007). Esta curva se compara luego con las 
curvas teóricas de los modelos de flujo pistón 
(REP), flujo pistón con dispersión (RFPD) o 
una combinación de ellos con reactores mezcla 
completa (RMC). 

Para la evaluación hidrodinámica se empleó 
la curva DTR siguiendo el método descrito en 
Levenspiel (2002). Como trazador se utilizó 
NaCl, inyectándose al sistema de manera 
instantánea a la entrada del humedal 300 g 
de la sal disueltos en dos litros de agua. Cada 
día se midieron la conductividad y caudal del 
efluente hasta observar los valores normales de 


(b) 


Figura 1. Fotografía del humedal construido con Cortadeira selloana: a) en 2010, b) en 2012. 
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conductividad del agua residual tratada. Para la 
medición de conductividad se utilizó una sonda 
multiparamétrica marca Horiba U-10. Luego se 
calculó el tiempo de retención hidráulico me- 
dio (t,,), según la ecuación (3), y se construyó 
la curva E normalizada (E,), normalizando el 
tiempo y la concentración del trazador mediante 
las ecuaciones (4) y (5): 


_ E(t,At,C,) 
" Y(at,C,) pen 


donde t, y C,son los tiempos y concentraciones 
de sal obtenidos en el transcurso del ensayo: 


(4) 


E, = 


C, 
SC, At, ha ds 


siendo 6 el tiempo normalizado. Por último se 
comparó la curva E, Us. O obtenida en el ensayo, 
con curvas E, normalizadas para patrones de 
flujo de REP ideal, REP con dispersión y con 
n reactores de mezcla completa conectados en 
serie o combinaciones de éstos. 


Evaluación cinética 


Los humedales artificiales pueden considerarse 
como reactores biológicos y su rendimiento 
ser estimado mediante una cinética de primer 
orden (ecuación (1)) para la remoción de DBO,, 
material en suspensión u otro contaminante, 
como el nitrógeno. Para la evaluación de la 
constante cinética se asumió una reacción de 
primer orden (Kadlec, 2003; Francia, Camarena, 
6 Guzmán, 2002; Reed et al., 1995) y se utilizaron 
las concentraciones de DBO, afluente y DBO, 
efluente obtenidos en ensayos realizados según 
APHA (APHA, AWWA, WEE, 2005), a partir de 
los 15 meses de construido el humedal, con 
una variación en el tiempo entre cada ensayo 
de unos tres meses (en septiembre y diciembre 


de 2011, y en marzo y mayo de 2012). Luego 
se validó dicha constante de reacción k con los 
datos registrados en 2010, cuando el caudal de 
alimentación era diferente (en 2010 era de 700 
1/d y desde febrero de 2011 fue de unos 222 1/d). 


Desarrollo de software de simulación 


Se desarrolló un programa de simulación escri- 
biendo un modelo matemático en el lenguaje 
Visual Basic 6.0 de Microsoft con las ecuaciones 
(1) y Q) (es decir, considerando que el sistema 
puede ser modelado como un RFP) y la cons- 
tante cinética obtenida. Este programa permite 
ensayar el comportamiento del sistema en estu- 
dio y de otros humedales análogos en diferentes 
situaciones. Permite simular escenarios para 
calcular y graficar la DBO, efluente, modifican- 
do cualquiera de las siguiente cinco variables: 
dimensiones L, W y h del humedal; porosidad 
e del medio; caudal diario del afluente v; DBO, 
afluente o las precipitaciones. Todas estas varia- 
bles afectan el tiempo de residencia del efluente 
en el sistema (8), y por lo tanto son determinan- 
tes en la DBO efluente y la capacidad de remo- 
ción del humedal. Para poder graficar la DBO, 
efluente en función de una de las variables, se 
deben introducir como datos las otras cuatro. 
Para simular la situación de lluvia se aplican 
las ecuaciones (6) y (7), a fin de contemplar el 
aumento de caudal y la dilución del afluente 
(considerando que toda el agua aportada por 
las lluvias ingresa a la entrada del humedal): 


A (LWh)e 
|, menlluavia LIV e 
1000 
DBO, afluente v 
DB Sy AAA 7 
O, afluente aparente mmlluvia LW ( ) 
1000 


donde t está expresado en días y DBO, afluente 
aparente es la DBO, que tendría el afluente si 
se diluyera con el agua aportada por la lluvia 
cada día. 
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Resultados y discusión 


Determinación de la 1 experimental del HAFS 
y modelo de flujo actual 


La figura 2 muestra la curva de distribución 
de tiempo de residencia realizada luego de 28 
meses de construido el humedal, en octubre de 
2012. A partir del análisis de la misma se obtuvo 
un f,, de 11.7 días. 

Comparando la curva E, obtenida (figura 2) y 
las curvas E, de reactores no ideales (Levenspiel, 
2002), se pudo ver la similitud de ésta con la cur- 
va E,, correspondiente a reactores con “flujo en 
caminos paralelos”. Advirtiendo esta similitud 
y observando en la curva E, obtenida la exis- 
tencia de dos picos acampanados, se propone 
un modelo hidrodinámico para el humedal en 
estudio, formado por dos RFP con dispersión en 
paralelo. El primer pico de la curva da indicios 
de un REP con flujo más rápido y, por su menor 
altura, una menor fracción de flujo respecto al 


segundo RFP más lento. Se propuso que cada 
uno de los RFP pueda ser representando por 
una función de distribución normal (f1(8) y 
f2(0)), caracterizada por sus parámetros media 
(que representarían el tiempo de residencia nor- 
malizado de cada reactor) y desviación estándar 
(que representaría la dispersión de cada uno de 
éstos). Se propuso que la suma ponderada de 
f1(0) y £2(0), según la fracción de flujo, dé como 
resultado una función (8) (ecuación (8)): 


a FLO)+b f2(0)= (0) (8) 


con a y b, factores de ponderación de las funcio- 
nes f1(0) y f2(0), respectivamente. Por último, 
se hace uso de la herramienta Solver de Excel 
para optimizar los parámetros de f1(0) y £2(0), 
a fin de obtener la mayor correlación (R?) entre 
la función (0) y la curva E, del humedal (figura 
3). Los parámetros de f1(0) y f2(0) buscados 
fueron los tiempos medios de residencia 0, y 0, 
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Figura 2. Curva E normalizada obtenida luego de 28 meses de construido el humedal. 
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(ecuaciones (9a) y (9b)), y desviaciones estándar 
O, Y 9, (ecuaciones (10a) y (10b)). 


9, - Ln (9a) 
En 
bra 
O, = 2 (9b) 
Eu 
D 
o, = [2 10a 
el ul ( ) 
D 
= 2 24 10b 
02 (10b) 


siendo D el coeficiente de dispersión longitudi- 
nal, u la velocidad media de flujo y L la longitud 
del reactor. 


Para el modelo mencionado de dos RFP con 
dispersión en paralelo se obtuvo un 6, de 0.550 
y un módulo de dispersión (D / ul) de 0.194 para 
el primer RFP y la fracción de flujo (parámetro a 
de la ecuación (8)) de 0.400; mientras que para el 
segundo REP, el 6, obtenido fue de 1.257; el mó- 
dulo de dispersión (D / ul) de 0.246 y la fracción 
de flujo (parámetro b de la ecuación (8)) de 0.600. 
El valor de R? encontrado en la optimación fue 
de 0.896, indicando una muy buena correlación 
entre la función (0) y la curva E, del humedal. 


Evaluación cinética. Determinación de la 
constante de reacción k para remoción de 
DBO, 


Con el tiempo de residencia medio experimental 
(11.7d) y la ecuación (1) se calculó la constante 
cinética, utilizando las concentraciones de DBO, 
afluente y DBO, efluente medidas entre sep- 
tiembre de 2011 y mayo de 2012. Con los valores 


el 


A 


s 


1 
! N 


82 


Figura 3. Modelo para evaluar la hidrodinámica del humedal; f1(8) y f2(0): funciones de distribución normal que representan 
los REP del modelo; 6, y 6,: tiempos de residencia normalizados para los REP. 
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de k calculados se determinó un valor promedio 
y se calculó la DBO, teórica del efluente, para 
determinar si este promedio puede representar 
el conjunto; caso contrario, se eliminan los va- 
lores con grandes desviaciones respecto al resto 
(cuadro 1). 

Se descartaron las determinaciones corres- 
pondientes al mes 15 (septiembre de 2011) por 
introducir una diferencia mayor a 25% entre 
la DBO, efluente medida y la estimada con el 
promedio de k. Por lo tanto, se determinó como 
constante de reacción k para la remoción de 
DBO, el promedio de las k calculadas para los 
meses 18, 21 y 23. Este valor resultó en k= 0.115 
1/d. 


Programa de simulación del HFSSH 


Las figuras 4 y 5 muestran la interfaz visual del 
software de simulación. Los rangos que admite 
el programa de simulación son los siguientes: 
un volumen de la celda (L* W* h) de 0 a 16 m?; 
una porosidad (e) de 0.3 a 0.6; un caudal diario 
(v) de 0 a 1 m*/d; DBO, afluente de 0 a 400 mg 
O,/1, y lluvias de 0 a 60 mm/d. Cabe destacar 
que el modelo sólo puede utilizarse en hume- 
dales que empleen para la remoción Cortaderia 
selloana, ya que ésta define la cinética del sistema 
y que el mismo no contempla variaciones de la 
remoción con la temperatura. Sin embargo, se 
podría hacer extensivo a otra vegetación (con 


Cuadro 1. Determinación de la constante k para remoción de DBO,. 


Edad del DBO, afluente | DBO, efluente k DBO, efluente teórica | DBO, efluente teórica 
humedal Fecha (mg/l) (mg/l) (1/d) | con k=0.105 1/d (mg/l) | con k=0.115 1/d (mg/D) 
15 meses Sep. 2011 162.8 66.6 47.5 Ñ 

18 meses Dic. 2011 147.5 39.93 43.0 38.4 

21 meses | Marzo 2012 92 24 26.8 24.0 

23 meses Mayo 2012 90.4 22.6 26.4 23.5 


k= 0.105 1/d, valor promedio de k considerando todas las determinaciones realizadas. 
k= 0.115 1/d, valor promedio de k eliminando las determinaciones realizadas en el mes 15 (septiembre de 2011). 


IM Simulador de HFSSH 


DBO Efluenze [m3 02/12] 


91,18 1 


Caudal Diano Mecio 102222 m3 / día 


Tiempo de Residencia Hidráulico 111,712 


DBO Afluente (123 mg 02/1 


Días .2,0€ 


72.9 


DBO Efluente 1313385  mg02/! 


Ancho 121 _M 
Arolo5 m| Vol Equivalente 2.,60ím3 pal 


Largo 15,1 m 


Porosidad [0,485 % 


2,2 


3,71 


43,59 


2€,47 


7,35 


Figura 4. Simulación de la variación de DBO, del efluente en función del caudal. 
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MÍ Simulador de HFSSH 


DSO Efluente (my 02/11 


22,59 1-- 
Caudal Diario Medio [0,2222 m3/día 


Tiempo de Residencia Hidráulico [11,712 — Días ES 
DBO Añuente [123 mg02/1 
DBO Efluente [31.385 mg02/1 


31,34 


esla. 
e 


Largo 15,1 m 
Porosidad |0.486 2% 


Precipitación [mm 110/dia) 


Figura 5. Simulación de la variación de DBO, del efluente en función de precipitaciones. 


otros valores de k) y contemplar la variación 
con la temperatura, introduciendo las variables 
correspondientes al programa. 


Evaluación de situaciones hipotéticas del 
HFSSH en estudio 


El límite máximo de DBO, efluente para la 
calidad de las aguas residuales para descarga 
a cursos de agua superficial, según las “Nor- 
mas para la protección de los recursos hídricos 
superficiales y subterráneos de la provincia de 
Córdoba — Decreto Provincial 415/99” es de 
30 mg O,/1. Para el sistema actual de volumen 
(LWh = 5.1 x 2.1 x 0.5) de 5.355 mi; e = 0.486; v 
= 0.22 m*/d; DBO, afluente = 123 mg O,/l y, 
en ausencia de lluvias, las simulaciones dan 
un DBO, efluente de 32 mg O,/], y no sería 
apto para descarga directa a un cuerpo de agua 
superficial. Se simularon situaciones variando 
cada una de estas variables para establecer los 
límites admisibles. 


Con precipitaciones 


La figura 6 muestra la variación de DBO, efluen- 
te simulado para diferentes situaciones de lluvia 


(mm de lluvia) para el humedal en estudio. El 
agua aportada por la lluvia tiene dos efectos 
encontrados en el sistema: a) un efecto desfa- 
vorable, pues disminuye el tiempo de residencia 
y así la remoción de materia orgánica, aumen- 
tando el DBO, efluente; y b) un efecto favorable, 
ya que disminuye la DBO, afluente aparente en 
situaciones de lluvia (ecuación (7)) por un efecto 
de dilución, disminuyendo la DBO, efluente. 
El primer efecto se observa hasta valores de 
precipitaciones de hasta 18 mm de lluvia; sin 
embargo, se observa que el efecto de la dilución 
prioriza a partir de ese valor, obteniéndose un 
efluente con una DBO, menor que los días sin 
lluvia. 


Con variación de DBO afluente 


La figura 7 muestra la simulación realizada, 
variando la DBO, afluente y manteniendo 
constante el resto de las variables (volumen de 
5.355 mí; e = 0.486; v = 0.22 m*/d) y ausencia 
de precipitaciones. El valor máximo de DBO, 
afluente admitido es de 115 mgO,/l, para cum- 
plir con el nivel de DBO, que exige el Decreto 
415-99 de 30 mg O,/l, para descarga a embalse 
(criterio para el cual fue diseñado el humedal). 
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Con variación del caudal para diferentes DBO, 
afluente 


La figura 8 ilustra la variación de la DBO, 
efluente en función del caudal diario v en caso 
de: a) mantenerse la DBO, afluente promedio; 
b) incrementarse en un 50%, y c) verse reducida 
un 50% para el humedal en estudio (volumen 
de 5.355 m?; e = 0.486) y en ausencia de precipi- 
taciones. En la figura 8 se observa que a medida 
que aumenta el caudal del afluente mantenien- 
do su DBO, constante, disminuye el tiempo 
de residencia y, por ende, baja la remoción de 
materia orgánica. Al variar la DBO, afluente, 
manteniendo el caudal constante, se ve afectada 
la DBO, efluente de igual manera, observándose 
que a mayores valores de DBO, afluente mayor 
es la DBO, efluente y cuanto menor es la DBO, 
en el afluente, menor en el efluente, como es de 
esperar. Se observa también que para cumplir 
con el nivel de DBO, efluente de 30 mgO, /1 que 
exige el Decreto 415-99 para una concentración 


de DBO afluente de 123 mgoO, / 1, el valor máxi- 
mo de caudal de entrada admitido es de 0.2123 
m?/d; para DBO afluente de 185.4 mgO,/1l es 
de 0.1647 m*/d, y para DBO afluente de 61.5 
mgoO,/l es de 0.4169 m*/d. 


Conclusiones 


Mediante la evaluación hidrodinámica se 
detectó que el sistema en estudio modificó 
el comportamiento de flujo pistón ideal y 
evolucionó hacia un modelo de dos reactores 
flujo pistón con dispersión en paralelo. Como 
se esperaba, el sistema ya en funcionamiento se 
alejó del modelo ideal, pues a través del tiempo 
se produce una variación del volumen efectivo 
del humedal y cambio en su morfología debido 
a poros ocupados por raíces de la vegetación 
y acumulación de sólidos, asentamiento de la 
grava y depresiones en el suelo por el peso de la 
celda, entre otros. Este sistema con canalizacio- 
nes, semejante a dos reactores de flujo pistón en 
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Figura 8. Variación de DBO, efluente en función del caudal V, para tres DBO, afluente. 
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paralelo con dispersión, se ajustó a un modelo 
Gaussiano. 

El cálculo de la constante cinética, tomando 
como modelo el RFP y las concentraciones de 
DBO, de la entrada y salida del humedal real, 
dio como resultado una k = 0.115 d”, aproxi- 
madamente siete veces menor a la constante 
utilizada para el diseño. Es decir, la eficiencia de 
remoción no coincidió con la supuesta al inicio 
para el diseño, tomada de la bibliografía. 

Estos dos factores, la modificación del com- 
portamiento hidrodinámico y el cálculo de la 
constante cinética real de cada humedal, resaltan 
la importancia de llevar a cabo estudios con 
sistemas a escalas reales y evaluarlos a través 
del tiempo. 

Si bien esta constante fue calculada median- 
te el modelo de flujo pistón ideal, es de gran 
utilidad para simular y evaluar los resultados 
en cuanto a la carga de materia orgánica en el 
efluente, factor primordial a controlar en efluen- 
tes cloacales que descargan a cuerpos de agua 
superficiales. La simulación permitió conocer 
a priori la concentración de materia orgánica 
del efluente del humedal para situaciones de 
interés, establecer los caudales y concentracio- 
nes de DBO, del afluente que puede tratar este 
humedal, así como el nivel de precipitaciones 
aceptables para el mismo, cumpliendo con los 
límites admisibles de la normativa local. 


Abreviaturas 


RFP: reactor flujo pistón. 


HES: humedal de flujo superficial. 

HFSS: humedal de flujo subsuperficial. 

HFSSH: humedal de flujo subsuperficial hori- 
zontal. 

HFSSV: humedal de flujo subsuperficial verti- 


cal. 

DBO, efluente: concentración de demanda 
bioquímica de oxígeno (DBO,) en el 
efluente. 

DBO, afluente: concentración de demanda 
bioquímica de oxígeno (DBO,) en el 
afluente. 

k constante de velocidad cinética. 


T tiempo de retención hidráulica. 

E: largo del humedal. 

W: ancho del humedal. 

h profundidad del humedal. 

€ porosidad del lecho. 

v: caudal volumétrico de efluente. 

DIPAS: Dirección Provincial de Agua y Sanea- 
miento. 

DTR: — Distribución de Tiempos de Residen- 
cia. 

LE tiempo de retención hidráulico medio. 

E,: curva E normalizada. 

t, tiempo l. 

C; concentración de sal obtenidos en el 
tiempo 1. 

6: tiempo normalizado. 

APHA: American Public Health Association. 


DBO, afluente aparente: DBO, que tendría el 
afluente si se diluyera con el agua 
aportada por la lluvia diariamente. 

fL(0), f2(0): funciones de distribución normal 

de los dos RFP con dispersión en 
paralelo. 

fracciones de flujo de los RFP 1 y 2 en 

paralelo, respectivamente. 

función f. 


0, 0, tiempos medios de residencia de los 
RFP 1 y 2 en paralelo, respectivamen- 
te. 

Or Op)! desviaciones estándar de los RFP 1 y 
2 en paralelo, respectivamente. 

D: coeficiente de dispersión longitudinal. 

u: velocidad media de flujo. 

D/ul: módulo de dispersión. 
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